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We generalize the recently proposed Greenberger-Horne-Zeilinger (GHZ) tripartite protocol [A. 
Galiautdinov, J. M. Martinis, Phys. Rev. A 78, 010305(R) (2008)] to fully connected networks of 
weakly coupled qubits interacting by way of anisotropic Heisenberg exchange g (XX + YY) +gZZ. 
Our model adopted here differs from the more familiar Ising-Heisenberg chain in that here every 
qubit interacts with every other qubit in the circuit. The assumption of identical couplings on 
all qubit pairs allows an elegant proof of the protocol for arbitrary N . In order to further make 
contact with experiment, we study fidelity degradation due to coupling imperfections by numerically 
simulating the N = 3 and N = 4 cases. Our simulations indicate that the best fidelity at unequal 
Oh couplings is achieved when (a) the system is initially prepared in the uniform superposition state 

(similarly to how it is done in the ideal case), and (b) the entangling time and the final rotations 
on each of the qubits are appropriately adjusted. 
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I. INTRODUCTION 



In the rapidly developing field of quantum computing the concept of quantum entanglement is considered to be 
of greatest practical importance. Many important applications,such as quantum communication secret sharing 
, open-destination teleportation [H , fault-tolerant computing 0, |j| , and others [(| , rely heavily on an architecture's 
ability to generate multipartite entangling states. This is especially true for the general Greenberger-Horne-Zeilinger 
(GHZ) states @, || 

m : 1 

|GHZ)^ = _(|0>®" + |1>®"), (1) 

CN ■ 

C\| ' _ 

also known as cat states (see, e.g., [9j), which, according to most entanglement measures, are maximally entangled 

' Various approaches to generating GHZ states have been proposed in the literature. A scheme f or g enerating cat 
states of a single- mode optical field by means of conditional measurement was proposed in Ref. [10(. A one-step 
multi-atomic GHZ state generation in a non-resonant cavity by way of cavity-assisted collisions was considered in Ref. 
PH . Ref. [IH discusses the possibility of generating GHZ states of electron spin qubits in a chain of quantum dots 
using the naturally available single-qubit rotations and the two-qubit Heisenberg exchange interaction. The minimum 
number of required operations in that proposal scales linearly with the number of qubits. In Ref. 13], in the context 
of capacitively coupled superconducting phase qubits, an Af-qubit GHZ state is proposed to be generated by applying 
an initial one-qubit Tr/2-pulse followed by an alternating sequence of 2(N — 2) + 1 controlled-NOT (CNOT) and SWAP 
gates. In Ref. [14| it is shown how the GHZ state can be generated in a multi-qubit system by coupling it sequentially 



to a mesoscopic static completely mixed spin bath. In Ref. 15|, in the context of circuit quantum electrodynamics, 
high-fidelity GHZ states (in a system whose mutual qubit detunings are large compared to the inter-qubit coupling 
strengths) are proposed to be probabilistically generated using one-qubit rotations and a single A^-qubit dispersive 
readout. 

In this paper we describe how GHZ states can be generated using a single ./V-qubit entangling pulse in a fully 
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connected network of qubits (see also [l6j]) interacting by way of anisotropic Heisenberg exchange [l7| 

N-l N 

^t = ^£ £ [9(44+<^) + ~gai^], (2) 
1=1 k=t+i 



where er^, \l = x,y,z, are the Pauli matrices ct m acting on the £-th qubit. We concentrate on the case g ^ g only, 
which in the context of Josephson phase qubits (our initial motivation [l8j]) corresponds to either capacitive (g = 0) 
or inductive (0 < \g/g\ < 0.1) coupling scheme. We show that the duration of the entangling pulse is given by 

teHZ = 7r/%-5l, (3) 

and thus in a typical experiment with g ~ 10 MHz the time iQHZ is on the order of 25 ns, which is small compared 
to the usually achieved coherence times of ~ 500 ns [lj| . We point out that the interaction Hamiltonian given in Eq. 
([2|) is the result of first projecting the exact system Hamiltonian onto the qubit subspace, switching to the interaction 
picture, and then applying the rotating wave approximation (RWA) in which the fast oscillating terms are ignored. 
Here, the time scale is set by a typical time to do a qubit operation, such as, e.g., the entangling time given above. 
The fast oscillations are the ones that occur at qubit transition frequencies, usually on the order of 10 GHz, and thus 
having transition times on the order of 0.1 ns [17] . 
Why employ fully connected networks? 

If we adopt an important simplifying assumption of identical couplings on all qubit pairs, then the answer to the 
above question is: "symmetry, simplicity, and experimental practicality." As will be shown below, the symmetric 
coupling scheme allows a rigorous proof of the protocol for arbitrary N. The resulting GHZ sequence is, indeed, very 
simple: prepare the system in the uniform superposition state, turn the coupling on and wait for a given time tcHZ) 
then apply a final, corrective iV-qubit rotation. 

There is a nice experimental realization of fully connected networks involving capacitively coupled superconducting 
qubits proposed by Matthew Neeley and his collaborators at UCSB [2(J. It consists of a superconducting island 
(an isolated piece of aluminum) connected via a capacitor to each of the qubits in the system. Using the standard 
methods of circuit analysis, one can transform this "star" configuration into an equivalent "delta" (or, fully connected 
network) configuration, in which there is a mutual capacitance between each pair of qubits. [The term "delta" comes 
from the three-qubit case, where the system looks like a triangle, a special case of the "star-delta" or "Y-dclta" 
transform in electrical engineering [ll[.] It is then straightforward to show that N qubits coupled to an island with 
capacitance C* are equivalent to N qubits coupled to each other in a complete graph with capacitance Ca = C+/N 
[20(. Thus, from the practical point of view, implementing a fully connected network using a solid state quantum 
computing architecture is relatively straightforward. There are still some residual errors in the coupling that are 
due to imperfections in the coupling capacitors, so the resulting pairwise couplings are not strictly identical. That 
somewhat limits applicability of the model and, in a rigorous analysis, must be taken into account. Below we show 
how the errors due to imperfections in the individual couplings can be drastically reduced by adjusting the entangling 
time and the final rotations in the N = 3, 4 cases. 

Our final motivation for the use of fully connected networks comes from their applicability to problems in number 
theory. Two of us (MWC and RD) have recently proposed an effective way to perform primality testing using this 
model with a modest number of qubits [22j • 

We now turn to the description and the proof of the GHZ protocol itself. 



II. GHZ PROTOCOL 



First, notice that for any number of qubits N the Hamiltonian @ has eigenenergy (see Appendix |A"| 

E Gnz = X (N) = C 2 N (g/2) (4) 

associated with the eigenstates |0)® and jl)®^, where Cjj — (N — l)N/2 is the total number of pairwise couplings 
present in the network. Then, the two families of the GHZ protocols are 

N = 2,4,6,... : e - iXotGSZ e^ N ^- x ^\GRZ) N = R 3 RiU ent Ri\0)® N , (5) 

and 

N = 3,5,7,... : e- l ^ tGHZ e l( -^ iN ^ )/2(7T /^\GRZ) N = R 2 U cnt R 1 \0) (SN , (6) 
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where 

R 3 = e-^+C-D^Vi/*, Rl = ®f =1 e— </ 4 , i? 2 = ®f =1 e"^/ 4 , (7) 

are the corresponding local rotations, and t/ ent = e~ iHinttGHZ is the entangling pulse of duration £ghz = 7r/2|<? — <jf|. 
To see in detail how these protocols work we first re-write the fully uniform superposition state, 

|^unifor m = J Rl|0) 8JV =^ 7I Yl l£l>"--l^>> (8) 

€ fc e{o,i} 

as a sum of partial uniform superpositions, each of which is characterized by the total number j of up spins in the 
respective direct-product components, as follows: 

1 N r— 

\^)uniior m =^J2V C N\ W ^ ( 9 ) 
3=0 

Here, \Wj) stands for a generalized VF-state (with j spins up), and C J N — N\/j\(N — j)\ is the corresponding binomial 
expansion coefficient. For example, for N = 3, \Wi) = [|001> + |010) + |100)]/\/3, \W 2 ) = [jOll) + |101) + |110)]/-y/3, 
etc. Notice that for any N, 

[\W ) + \W N )]/V2 = [\0)® N + \1}® N ]/V2 = \GRZ) N . (10) 

We now state an important property of the interaction Hamiltonian -ff; n t: 
The eigenvalue formula: For any j = 0, 1, 2, . . . , N, the states \Wj) are the eigenstates of H lnt associated with the 
eigenvalue 

Xj(N) = j(N - j)(g - g) + C 2 N {g/2). (11) 

Proof outline. The operator S z = J^ fe a k z commutes with Hi nt , providing a good quantum number, the total number 
of up spins in a given state. Since i?; n t and T, z share eigenvectors, the \Wj) are eigensubspaces. The proof of the 
analytical formula (TTTj) for the eigenvalues is given in Appendix [A] Notice that since the Hamiltonian |2]) commutes 
with the operator X® N = NOT®^, if \Wj) is an eigenvector, then so is \Wn-j), with the same eigenvalue. This is a 
powerful property, as it immediately tells us that the eigenenergies are constrained to satisfy the relation Xj = Ajv-j- 
■ 

The eigenvalue formula immediately allows us to write down the effect of the entangling pulse on the uniform 
superposition, 

i N r— 

t/ent ^uniform = ^ E V °* e^^Wj) , (12) 
3=0 

giving the state on which the final A-qubit rotations R3R1, R2 will be acting. Assuming g > g, we have the useful 
relation, 

g-iVoHz = e - iX ° tGS * exp{-i\j(N - j)(g - g)t G Hz]} = e - iX ° tasz exp[-i(ir/2)j(N - j)} = e - iX ° tGHZ (~iy (N - j \ (13) 
Now, for N odd (Appendix |B|) . 

i? 2 - 1 |GHZ) N odd = ^72 E + y )HN] \ W i)- ( 14 ) 

3=0 v 

In this case i N — ±i. Recalling that 

(W j \W k )=6 jk , (15) 
we see that the full state is properly normalized. Combining Eqs. (TT2")) . (TT3"|) , p^|) . with the identity 

iW-fl = (-i)H N -n = cos\j(N -j)ix/2] 

= ex P [a-ir^] ^[1 + i-iyi»] = ex P ^±^m , ( 16) 
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gives the GHZ state protocol for odd N. 
For N even (Appendix 

i N I — -f_i\j I p -iB 

i?r 1 [i? 3 ]^ 1 |GHZ) N cvcn = —^0%^-^ \ Wj ), (17) 

3=0 V 

where 

6(N) = (*/2)[2+(-l) lf ' 2 ]. (18) 
Combining Eqs. (fT2|) and (fT7|) with the identity 

(-iy(N-3l = e *»Wa-i)£^![(_iy + e -*], (19) 
v2 

gives the GHZ state protocol for even A. [When g < g, and TV is even, an extra factor of e - i(7r/4)(<T ' +a " ) should be 
inserted after i?3, otherwise the GHZ state accumulates an internal phase of (-1).] 

Notice that for g = g, the uniform superposition state is an eigenstate of the interaction Hamiltonian Hi nt . Thus, 
the | V') uniform does not change (apart from accumulating an overall phase) under the action of the entangling pulse 
U cn t- This explains the previously stated requirement g ^ g (cf. [13 )• 



III. FIDELITY DEGRADATION DUE TO COUPLING ERRORS 



Achieving identical couplings on all qubit pairs is experimentally difficult. One way to reduce errors due to coupling 
imperfections is to adjust the final qubit rotations and the entangling time in such a way as to maximize the fidelity 
of the resulting state. Figure [1] shows the result of numerical optimization with respect to the Frobenius distance 
between the generated state and the GHZ target, 



T=l- J Mo* - (GHZ|)(|V>o P t - |GHZ», 



(20) 



for A~ = 3. Here the Hamiltonian is given by 

(.912/2) [{alal + affi) + (1 - ry 23 ) (a 2 x a\ 



+ 0*0*) + (1 - 7] 13 ) {a l x al 



^-2^-3 

a, a. 



alal) 



where 1712 > is the reference coupling on qubits 1 and 2, 7723, 7713 > are the coupling errors on pairs 23 and 13, and n 
is the ZZ coupling, which for simplicity is assumed to be constant throughout the network (the results of this section 
do not depend on this assumption). Our simulations indicate that the best fidelity is achieved when the initial pulse 
is chosen to generate the uniform superposition state, similar to how it is done in the ideal case with 7723 = 7713 = 0. 
A typical optimized state, here shown at k = 0.05, 7723 = 0.02, and 7713 = 0.06, has the form 

/ 0.707099 \ 
0.000692 
-0.001956 
0.002566 
/ " pl ~ 0.002566 
-0.001956 
0.000692 
V 0.707099 / 



(22) 



which corresponds to i pt/(7r/2(7i 2 (l — n)) 
be compared to the non-corrected state, 



1.0505, a^'VU/ 2 ) = 0.9785,0.9713,0.9825, and T = 0.9953. This may 



/ 0.706616 \ 

0.000697 +0.01543H 



/ original 



-0.002792 
0.002988 4 
0.002988 4 

-0.002792 



- 0.010929?: 
0.017844i 
0.017844i 

- 0.010929?: 



0.000697 4-0.015431?; 
0.706616 



(23) 
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FIG. 1: (Color online) Optimal values of the final rotation angles a* pt , i — 1,2, 3, and the entangling time t op t that minimize 
the error due to coupling imperfections in the N = 3 case. The gate fidelity is defined in Eq. f|20 J> . The corresponding 
Hamiltonian is given in Eq. (|21[) . Here, k — 0.05, 7723 = 0.02, and < 7713 < 0.10. Fidelity curve for the original protocol run 
with no corrections is also shown. 



that has T — 0.9628, which we get by running the original protocol at the same values of k and 7723, 7713. 

We have also performed optimization in the N = 4 case (not shown) where it was found that in order to maximize 
the fidelity of the resulting state only the entangling time tcnz and the final Z-rotation on the first qubit have to be 
adjusted. The R± pulse appearing on both sides of U ea t should, however, remain the same. 

The following argument indicates that similar optimization results may be expected for arbitrary N. Assuming 
coupling errors in the interaction Hamiltonian to be sufficiently small, their effect is to perturb (in general all of) the 
N + 1 eigenvalues Xj corresponding to the \Wj) states. If we then optimize the fidelity on the evolution time £ghz 
and iV final single-qubit rotation angles, we have N + 1 parameters to do so. We thus have N + 1 parameters in 
which to correct N + 1 state coefficient errors due to N + 1 perturbed phases coming from the perturbed eigenvalues. 
Minimizing the Frobenius distance needed for the fidelity is essentially solving a nonlinear least squares problem. 
Under only weak conditions, we will have a solution guaranteed as we have N + 1 parameters at our disposal. If 
there were errors in the initial fully uniform superposition state, then the N + 1 parameters would be insufficient in 
number to correct errors from both the initial state and a perturbed Hamiltonian. Improvements would be possible, 
but the correct initial state would then seem to be the best place to start from when precisely N + 1 parameters are 
available to later correct errors in the protocol. This argument shows that finding corrected fidelity by adjusting the 
entangling time and the final N rotation angles is always possible. It does not state how well the correction can be 
done, but the N = 3 and N = 4 cases presented above give illustrations. 
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IV. CONCLUSION 



In summary, we have demonstrated how the single-step GHZ protocol proposed for triangularly coupled supercon- 
ducting qubit system can be extended to a fully connected network of an arbitrary number of qubits. Two types of 
protocols can be identified here, depending on whether N is odd or even. When N is odd, the protocol (up to an 
overall insignificant phase) is the same as in the three-qubit case; it is given by a generic sequence A^t/ent^^lO)®^, 
with and representing the corresponding x and y rotations performed on the N qubits. When N is even, 

the final rotation gets modified and the sequence becomes ^5 / 7 r/2f4nt5 / 7r/2|0)® Ar , where Z\ is a z-rotation on the first 
qubit by angle which is given in Eq. (|18p . As a subproblem to describe the dynamics in the \Wj) subspaces, we 
exactly solved for the corresponding eigenenergies. 

We emphasize again that implementing identical couplings on all pairs of qubits is a difficult task. Here we 
proposed a possible approach to reducing the errors due to coupling imperfections in the experimentally important 
low-dimensional cases. 



APPENDIX A: PROOF OF THE EIGENVALUE FORMULA, EQ. (fTTj) . 



1. Raising and lowering operators 



For the purposes of this Appendix the matrix representations of states |0) = (0, 1) T and |1) = (1, 0) T are taken. 
We first define the raising and lowering operators S± = | (<r* ± ia k ) with the following properties: 

S-\0) = 0, S+|l) = 0, (S-) 2 = (S+) 2 = 0, S+|0) = |1>, S-\l) = |0>, 

k ok ok k ok ok ok k ok ok k ok 

O ~ O I — O | . u ^ iD — O ■ 0_|_fT^ — lj i ? O (T^, - — O j 

a\ = h 2S k _S k + , a z \0) = -|0>, a z \l) = |1), 

and 



S^ 1 S^ 2 



It will be useful to define the operators 



N 



N 



k=l 



k=l 



obeying 



[£+,E_]=£ z , [£ z ,£±]=±2E d 



Using these commutation relations we arrive at the following useful products. 
Lemma 1: 

E,(E ± )'' = (£ ± )'£,±2j(£±)'\ 

Proof (by induction). For j = 1, 

[£z, E ± ] = E Z E± — £±£ z = ±2£±, E Z E± = £±£ z ± 2£±. 

Assuming that £ z (£±)" = (E±)"E Z ± 2n(E ± ) n , we have 

E.(E ± )" +1 = (E 2 (E ± )")E± 

= ((£±) n £ z ± 2rc(£±) n ) £± 

= (E±) n (E±E, ± 2E±) ± 2n(£±)"+ 1 

= (£ ± )" +1 £ z ±2(n + l)(£ ± )" +1 . ■ 

Lemma 2: 



(Al) 
(A2) 



(A3) 



(A4) 



(A5) 



E_(E+y = (E+yS_ - jCE+y- 1 ^ - j(j - 1)(E+) 



(A6) 
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Proof (by induction). For j = 1, 

For j = 2, 

= £ + (£ + £_ — Tj z ) — — 2E + 

= (S+) 2 E_ - 2S+S Z - 2S+. 

Now, assuming £_(£+)" = (£+)"£- - - n(n - we get 

£_(£ + )" +1 = (£_(£+)")£+ 

= (£+)"(£_£+) - - n(„ - 1)(£ + )" 

- (£+)" +1 £_ - (£+)"£, - n(£+)"£ z - 2n(£+)" - n(n - 1)(£+)" 

= (£+)" +1 £_ - (n + !)(£+)"£* - (n(n - 1) + 2n)(S+) n 

= (£ + )" +1 £_ - (n + l)(S+rE, - ((n + l)((n + 1) - l))(E+) n . ■ 

We now generate partial uniform superpositions (generalized W states) using the raising and lowering operators: 

I Wb) = 1|0)® W , 



\Wi) 

\w 2 ) 



1 N 

=E^io>^, 



'C 



N fe=l 

JV-1 W 



y iV fe 1= l fe 2 >fci 



it^) = ^L= E E - E E si^-st'-^w 



C^j- fci = l k 2 >k\ kj-\>kj-2 kj>kj-\ 

If wc take into account all double counting we can see that these states can also be written in a form 
| Wb) = l|0) 8Ar , 



\ Wl ) = -±= £S*|0)®" = -1=£ + |0>^, 

V °JV fe=l V C 7V 

|W 2 ) 



V AT fc 1= lfe 2 =l 



(A7) 



Thus, 



N N N N 



" ^=EE _ _ 

p-yC 3 N fei=i fe 2 =i /£,_!=! fc,=i 



/ S^ 1 5^...S+ i - 1 5+ < |0)® JV = =(£+^10)®^. 



(A8) 



S + |^) 



(£ + ) 3+i |0)® iV , 



-^=fo' + i)i>/c^ r )|Wi + i). 



= A/(JV-j)(j + l)|W,-+i). 



(A9) 



Using the lemmas and the facts that E.lO)®^ = and E^O)®^ = -N\0)® N , we find: 
£_|I^} = — ^(ME+niO)^, 



"'"V+r'io) 



J-\ 



and 



—.((n + yx z + 2j(z + y)\o)® N , 



(-N + 2j)(X + y\0)® N , 



3'y 
1 



= (2j-N)\W j ). 

2. Eigenvalues of H lnt 



First notice that 



from which it follows that 



^Nx 2 + Vy'Vy 2 = 2 (S^S* 2 + St 1 Sty , 

H -t = ^E E (^(s^S^+S^sty+g^aty 



ki —1 k^ — k\ +1 

Also, 

ki,k2 — l ki^k2 — l k—1 fci=l fc 2 =/ci+l 

and similarly, 

TV N N N-l N N 

= s + s - 2 = E s + 15 '- 2 + E 5 + 5 '- = 2 E E s + s - 2 + E 5 + 5 '- 
fei,fe 2 =i fci^fe 2 =i fe=i fci=ife 2 =fei+i fe=i 

TV AT AT JV-1 N N 

= s - s + = E s '- 15 '+ + E s - s + = 2 E E s '- 15 + 2 + E s - s + 
fci,fe 2 =i *;i^*; 2 =i fe=i fei=ife 2 =fei+i fe=i 

Using the fact that J2k=i S + S - + EfeLi 5 '-' S '+ = ^ we re-write # int in the form 

/ N-l N N-l N N-l N \ 

^= %E E s^+asE E st^+ffE E <^ fe2 

V fei=ife 2 =fei+i fei=ife 2 =fei+i fci=ife 2 =fci+i / 

= 1 (V 9 ^(E+S- - E + 2 -9 - E + .9 - NI) 

= 1 (2 5 (E+£_ + + . 9 S Z S, - (2 5 + g)7VI) , 
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tf int = i(2 ff tff nt +p4), (A14) 

where 

Hf nt = + £_E+ - NI, Hf nt = S Z S Z - JVI. (A15) 

We now act with f/^ t on the partial uniform states, 

HPjWj) = (£+£- + — NT) \Wj), 



= VJ(N - J + + V(N - j)(J + m-\W j+1 ) - N\Wj), 

= U(N - j + 1) + (N - j)(j + 1) - N) \W 3 ), 

= 2j(N-j)\Wj), (A16) 

and similarly, 

HfjW,) = - NT)\W 3 ) = {{2j - Nf - N)\W 3 ) (A17) 

It is then easy to see that H^Wj) = \j(N)\Wj), where 

Xj(N) = \ (4j(N - j)g + ((2j - N) 2 - N)g) = j(N - j)(g ~ g) + C 2 N (g/2) . (A18) 
This completes the proof of the eigenvalue formula. 

APPENDIX B: PROOF OF EQ. CHD 

Here, N is odd. Using the fact that for y S {0, 1}, 

-i7rtr a ,/4| 



|y) = ^[|0)- i (-lf|l>] = ^ J2 H)*(-l) w 1*>, (Bl) 
we see that the action of R^ 1 on arbitrary direct-product states is given by 

jVl+V2 + ---+VN 

R^\yiy 2 ...yN) = ^ E z Zl+ ^ + - +2 "(-l)^ Zl+ ^ + -+^ z ~|z 1 )|z 2 ) . . . |**>. (B2) 

z k e{0A} 

We now let Ztot = J3fcLi 2 fc when occurring in such sums, and write more compactly 

jVl+V2 + ---+VN 

R^\y lV2 ...y N ) = J2 i Ztat (-inM...\^N). (B3) 

2fcG{0,l} 

We note that < Ztot < We then have 

1 ,-JV 

A^O)®* - ^72 E ^ l*i> ■ • ■ M' ^11)^ = E ^(-1)^1^1) • • ■ M- (B4) 

Therefore, 

i? 2 - 1 |GHZ) N odd = i?- 1 -L(|0)^ + |1>®") = -j^j- 2 E^[l + (-1)^]^) . . . |^). (B5) 

Z 3 

Replacing the sum over z^s by z tot , which is just an integer j in the range 0, . . . , N, gives 

i? 2 - 1 |GHZ) N odd = ^Ev / ^ P[1+ y )J ' Ar V 3 )- (B6) 
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APPENDIX C: PROOF OF EQ. (fTf|) 

Here, N is even. Note that for x € {0, 1}, 



which gives 



Since 



e --v/4| x)= i [( _ 1)10) + | 1)]= i j2 (-ir^ 

V2 V2 ,e{o,i} 



R i 1 \ x ) = ^Wj2 E (-i)^(-i)^ki) 



z 2 ) • ■ ■ \z N ) 



ZjS{0,l} 



(CI) 



(C2) 



(C3) 



we have 



J Rr 1 [^ 3 ]- 1 |GHZ) Ncvon = i?r 1 ^|(|0)® Ar + e^|l)® JV ) 



[(-lJJ+e**] 



where y> = —6(N). 



j=Q V 
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